Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs.

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

Base editing strategies to convert CAG to CAA diminish the disease-causing mutation in Huntington's disease.

An expanded CAG repeat in the huntingtin gene ( HTT ) causes Huntington's disease (HD). Since the length of uninterrupted CAG repeat, not polyglutamine, determines the age-at-onset in HD, base editing strategies to convert CAG to CAA are anticipated to delay onset by shortening the uninterrupted CAG repeat. Here, we developed base editing strategies to convert CAG in the repeat to CAA and determined their molecular outcomes and effects on relevant disease phenotypes. Base editing strategies employing combinations of cytosine base editors and gRNAs efficiently converted CAG to CAA at various sites in the CAG repeat without generating significant indels, off-target edits, or transcriptome alterations, demonstrating their feasibility and specificity. Candidate BE strategies converted CAG to CAA on both expanded and non-expanded CAG repeats without altering HTT mRNA and protein levels. In addition, somatic CAG repeat expansion, which is the major disease driver in HD, was significantly decreased by a candidate BE strategy treatment in HD knock-in mice carrying canonical CAG repeats. Notably, CAG repeat expansion was abolished entirely in HD knock-in mice carrying CAA-interrupted repeats, supporting the therapeutic potential of CAG-to-CAA conversion base editing strategies in HD and potentially other repeat expansion disorders.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements.

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Huntingtin turnover: modulation of huntingtin degradation by cAMP-dependent protein kinase A (PKA) phosphorylation of C-HEAT domain Ser2550.

Huntington's disease (HD) is a neurodegenerative disorder caused by an inherited unstable HTT CAG repeat that expands further, thereby eliciting a disease process that may be initiated by polyglutamine-expanded huntingtin or a short polyglutamine-product. Phosphorylation of selected candidate residues is reported to mediate polyglutamine-fragment degradation and toxicity. Here to support the discovery of phosphosites involved in the life-cycle of (full-length) huntingtin, we employed mass spectrometry-based phosphoproteomics to systematically identify sites in purified huntingtin and in the endogenous protein by proteomic and phosphoproteomic analyses of members of an HD neuronal progenitor cell panel. Our results bring total huntingtin phosphosites to 95, with more located in the N-HEAT domain relative to numbers in the Bridge and C-HEAT domains. Moreover, phosphorylation of C-HEAT Ser2550 by cAMP-dependent protein kinase (PKA), the top hit in kinase activity screens, was found to hasten huntingtin degradation, such that levels of the catalytic subunit (PRKACA) were inversely related to huntingtin levels. Taken together, these findings highlight categories of phosphosites that merit further study and provide a phosphosite kinase pair (pSer2550-PKA) with which to investigate the biological processes that regulate huntingtin degradation and thereby influence the steady state levels of huntingtin in HD cells.

Mutations causing Lopes-Maciel-Rodan syndrome are huntingtin hypomorphs.

Huntington's disease pathogenesis involves a genetic gain-of-function toxicity mechanism triggered by the expanded HTT CAG repeat. Current therapeutic efforts aim to suppress expression of total or mutant huntingtin, though the relationship of huntingtin's normal activities to the gain-of-function mechanism and what the effects of huntingtin-lowering might be are unclear. Here, we have re-investigated a rare family segregating two presumed HTT loss-of-function (LoF) variants associated with the developmental disorder, Lopes-Maciel-Rodan syndrome (LOMARS), using whole-genome sequencing of DNA from cell lines, in conjunction with analysis of mRNA and protein expression. Our findings correct the muddled annotation of these HTT variants, reaffirm they are the genetic cause of the LOMARS phenotype and demonstrate that each variant is a huntingtin hypomorphic mutation. The NM_002111.8: c.4469+1G>A splice donor variant results in aberrant (exon 34) splicing and severely reduced mRNA, whereas, surprisingly, the NM_002111.8: c.8157T>A NP_002102.4: Phe2719Leu missense variant results in abnormally rapid turnover of the Leu2719 huntingtin protein. Thus, although rare and subject to an as yet unknown LoF intolerance at the population level, bona fide HTT LoF variants can be transmitted by normal individuals leading to severe consequences in compound heterozygotes due to huntingtin deficiency.

Promotion of somatic CAG repeat expansion by Fan1 knock-out in Huntington's disease knock-in mice is blocked by Mlh1 knock-out.

Recent genome-wide association studies of age-at-onset in Huntington's disease (HD) point to distinct modes of potential disease modification: altering the rate of somatic expansion of the HTT CAG repeat or altering the resulting CAG threshold length-triggered toxicity process. Here, we evaluated the mouse orthologs of two HD age-at-onset modifier genes, FAN1 and RRM2B, for an influence on somatic instability of the expanded CAG repeat in Htt CAG knock-in mice. Fan1 knock-out increased somatic expansion of Htt CAG repeats, in the juvenile- and the adult-onset HD ranges, whereas knock-out of Rrm2b did not greatly alter somatic Htt CAG repeat instability. Simultaneous knock-out of Mlh1, the ortholog of a third HD age-at-onset modifier gene (MLH1), which suppresses somatic expansion of the Htt knock-in CAG repeat, blocked the Fan1 knock-out-induced acceleration of somatic CAG expansion. This genetic interaction indicates that functional MLH1 is required for the CAG repeat destabilizing effect of FAN1 loss. Thus, in HD, it is uncertain whether the RRM2B modifier effect on timing of onset may be due to a DNA instability mechanism. In contrast, the FAN1 modifier effects reveal that functional FAN1 acts to suppress somatic CAG repeat expansion, likely in genetic interaction with other DNA instability modifiers whose combined effects can hasten or delay onset and other CAG repeat length-driven phenotypes.

Genetic and Functional Analyses Point to FAN1 as the Source of Multiple Huntington Disease Modifier Effects.

A recent genome-wide association study of Huntington disease (HD) implicated genes involved in DNA maintenance processes as modifiers of onset, including multiple genome-wide significant signals in a chr15 region containing the DNA repair gene Fanconi-Associated Nuclease 1 (FAN1). Here, we have carried out detailed genetic, molecular, and cellular investigation of the modifiers at this locus. We find that missense changes within or near the DNA-binding domain (p.Arg507His and p.Arg377Trp) reduce FAN1's DNA-binding activity and its capacity to rescue mitomycin C-induced cytotoxicity, accounting for two infrequent onset-hastening modifier signals. We also idenified a third onset-hastening modifier signal whose mechanism of action remains uncertain but does not involve an amino acid change in FAN1. We present additional evidence that a frequent onset-delaying modifier signal does not alter FAN1 coding sequence but is associated with increased FAN1 mRNA expression in the cerebral cortex. Consistent with these findings and other cellular overexpression and/or suppression studies, knockout of FAN1 increased CAG repeat expansion in HD-induced pluripotent stem cells. Together, these studies support the process of somatic CAG repeat expansion as a therapeutic target in HD, and they clearly indicate that multiple genetic variations act by different means through FAN1 to influence HD onset in a manner that is largely additive, except in the rare circumstance that two onset-hastening alleles are present. Thus, an individual's particular combination of FAN1 haplotypes may influence their suitability for HD clinical trials, particularly if the therapeutic agent aims to reduce CAG repeat instability.

Genetic Risk Underlying Psychiatric and Cognitive Symptoms in Huntington's Disease.

BACKGROUND: Huntington's disease (HD) is an inherited neurodegenerative disorder caused by an expanded CAG repeat in the HTT gene. It is diagnosed following a standardized examination of motor control and often presents with cognitive decline and psychiatric symptoms. Recent studies have detected genetic loci modifying the age at onset of motor symptoms in HD, but genetic factors influencing cognitive and psychiatric presentations are unknown. METHODS: We tested the hypothesis that psychiatric and cognitive symptoms in HD are influenced by the same common genetic variation as in the general population by 1) constructing polygenic risk scores from large genome-wide association studies of psychiatric and neurodegenerative disorders and of intelligence and 2) testing for correlation with the presence of psychiatric and cognitive symptoms in a large sample (n = 5160) of patients with HD. RESULTS: Polygenic risk score for major depression was associated specifically with increased risk of depression in HD, as was schizophrenia risk score with psychosis and irritability. Cognitive impairment and apathy were associated with reduced polygenic risk score for intelligence. CONCLUSIONS: Polygenic risk scores for psychiatric disorders, particularly depression and schizophrenia, are associated with increased risk of the corresponding psychiatric symptoms in HD, suggesting a common genetic liability. However, the genetic liability to cognitive impairment and apathy appears to be distinct from other psychiatric symptoms in HD. No associations were observed between HD symptoms and risk scores for other neurodegenerative disorders. These data provide a rationale for treatments effective in depression and schizophrenia to be used to treat depression and psychotic symptoms in HD.

Acquisition of an oncogenic fusion protein is sufficient to globally alter the landscape of miRNA expression to inhibit myogenic differentiation.

The differentiation status of tumors is used as a prognostic indicator, with tumors comprised of less differentiated cells exhibiting higher levels of aggressiveness that correlate with a poor prognosis. Although oncogenes contribute to blocking differentiation, it is not clear how they globally alter miRNA expression during differentiation to achieve this result. The pediatric sarcoma Alveolar Rhabdomyosarcoma, which is primarily characterized by the expression of the PAX3-FOXO1 oncogenic fusion protein, consists of undifferentiated muscle cells. However, it is unclear what role PAX3-FOXO1 plays in promoting the undifferentiated state. We demonstrate that expression of PAX3-FOXO1 globally alters the expression of over 80 individual miRNA during early myogenic differentiation, resulting in three primary effects: 1) inhibition of the expression of 51 miRNA essential for promoting myogenesis, 2) promoting the aberrant expression of 43 miRNA not normally expressed during myogenesis, and 3) altering the expression pattern of 39 additional miRNA. Combined, these changes are predicted to have an overall negative effect on myogenic differentiation. This is one of the first studies describing how an oncogene globally alters miRNA expression to block differentiation and has clinical implications for the development of much needed multi-faceted tumor-specific therapeutic regimens.

Acquisition of an oncogenic fusion protein serves as an initial driving mutation by inducing aneuploidy and overriding proliferative defects.

While many solid tumors are defined by the presence of a particular oncogene, the role that this oncogene plays in driving transformation through the acquisition of aneuploidy and overcoming growth arrest are often not known. Further, although aneuploidy is present in many solid tumors, it is not clear whether it is the cause or effect of malignant transformation. The childhood sarcoma, Alveolar Rhabdomyosarcoma (ARMS), is primarily defined by the t(2;13)(q35;q14) translocation, creating the PAX3-FOXO1 fusion protein. It is unclear what role PAX3-FOXO1 plays in the initial stages of tumor development through the acquisition and persistence of aneuploidy. In this study we demonstrate that PAX3-FOXO1 serves as a driver mutation to initiate a cascade of mRNA and miRNA changes that ultimately reprogram proliferating myoblasts to induce the formation of ARMS. We present evidence that cells containing PAX3-FOXO1 have changes in the expression of mRNA and miRNA essential for maintaining proper chromosome number and structure thereby promoting aneuploidy. Further, we demonstrate that the presence of PAX3-FOXO1 alters the expression of growth factor related mRNA and miRNA, thereby overriding aneuploid-dependent growth arrest. Finally, we present evidence that phosphorylation of PAX3-FOXO1 contributes to these changes. This is one of the first studies describing how an oncogene and post-translational modifications drive the development of a tumor through the acquisition and persistence of aneuploidy. This mechanism has implications for other solid tumors where large-scale genomics studies may elucidate how global alterations contribute to tumor phenotypes allowing the development of much needed multi-faceted tumor-specific therapeutic regimens.

Familial co-segregation of Coffin-Lowry syndrome inherited from the mother and autosomal dominant Waardenburg type IV syndrome due to deletion of EDNRB inherited from the father.

We report an African-American family that was identified after the proposita was referred for diagnostic evaluation at 4½ months with a history of Hirschsprung and dysmorphic features typical of Waardenburg syndrome (WS). Family evaluation revealed that the father had heterochromidia irides and hypertelorism supporting the clinical diagnosis of WS; however, examination of the mother revealed characteristic facial and digital features of Coffin-Lowry syndrome (CLS). Molecular testing of the mother identified a novel 2 bp deletion (c.865_866delCA) in codon 289 of RPS6KA3 leading to a frame-shift and premature termination of translation 5 codons downstream (NM_004586.2:p.Gln289ValfsX5). This deletion also was identified in the proposita and her three sisters with a clinical suspicion of CLS, all of whom as carriers for this X-linked disorder had very subtle manifestations. The molecular confirmation of WS type 4 (Shah-Waardenburg; WS4) was not as straightforward. To evaluate WS types 1-4, multiple sequential molecular tests were requested, including Sanger sequencing of all exons, and deletion/duplication analysis using MLPA for PAX3, MITF, SOX10, EDN3 and EDNRB. Although sequencing did not identify any disease causing variants, MLPA identified a heterozygous deletion of the entire EDNRB in the father. This deletion was also found in the proposita and the oldest child. Since the heterozygous deletion was the only change identified in EDNRB, this family represents one of the few cases of an autosomal dominant inheritance of WS4 involving the endothelin pathway. Altogether, clinical evaluation of the family revealed one child to be positive for WS4 and two positive for CLS, while two children were positive for both diseases simultaneously (including the proposita) while another pair test negative for either disease. This kinship is an example of the coincidence of two conditions co-segregating in one family, with variable phenotypes requiring molecular testing to confirm the clinical diagnoses.

Identification of CK2 as the kinase that phosphorylates Pax3 at Ser209 in early myogenic differentiation.

The myogenic transcription factor Pax3, a member of the paired class homeodomain family of transcription factors, plays an essential role in early skeletal muscle development. We previously demonstrated that Pax3 is phosphorylated at three specific residues (Ser201, Ser205, and Ser209) and that the pattern of phosphorylation at these sites changes throughout early myogenesis. Further, we demonstrated that the protein kinase CK2 phosphorylates Pax3 at Ser205 and that this phosphorylation event is required for the subsequent phosphorylation of Ser201 by GSK3ß. However, the kinase that phosphorylates Pax3 at Ser209 has yet to be identified. In the present work we use standard purification methods and in vitro biochemical analyses to provide solid evidence identifying the protein kinase CK2 as phosphorylating Pax3 at Ser209. Further, we qualitatively demonstrate that the phosphorylation of Pax3 at Ser209 by CK2 is enhanced when Ser205 is previously phosphorylated. Taken together, our results allow us to propose a mechanism to describe the ordered phosphorylation of Pax3 throughout early myogenesis.

Identification of serines 201 and 209 as sites of Pax3 phosphorylation and the altered phosphorylation status of Pax3-FOXO1 during early myogenic differentiation.

Pax3, a member of the paired class homeodomain family of transcription factors, is essential for early skeletal muscle development and is key in the development of the childhood solid muscle tumor alveolar rhabdomyosarcoma (ARMS). ARMS is primarily characterized by a t(2;13)(q35;q14) chromosomal translocation, which fuses the 5'-coding sequences of Pax3 with the 3'-coding sequence of the forkhead transcription factor FOXO1 generating the oncogenic fusion protein Pax3-FOXO1. We previously demonstrated that Pax3 and Pax3-FOXO1 are phosphorylated by the protein kinase CK2 at serine 205 in proliferating primary myoblasts and that this phosphorylation event is rapidly lost from Pax3, but not Pax3-FOXO1 upon the induction of differentiation. However, reports suggested that additional sites of phosphorylation might be present on Pax3. In this report we use in vitro and in vivo analyses to identify serines 201 and 209 as additional sites of phosphorylation and along with serine 205 are the only sites of phosphorylation on Pax3. We provide solid evidence supporting the role of the protein kinase GSK3ß as phosphorylating Pax3 at serine 201. Using phospho-specific antibodies we demonstrate a changing pattern of phosphorylation at serines 201, 205, and 209 throughout early myogenic differentiation and that this pattern of phosphorylation is different for Pax3-FOXO1 in primary myoblasts and in several ARMS cell lines. Taken together, our results allow us to propose a molecular model to describe the changing pattern of phosphorylation for Pax3 and the altered phosphorylation for Pax3-FOXO1 during early myogenic differentiation.